Homogeneous MERA states: an information theoretical analysis 
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Homogeneous Multi-scale Entanglement Renormalization Ansazt (MERA) state have been re- 
cently introduced to describe quantum critical systems. Here we present an extensive analysis of 
the properties of such states by clarifying the definition of their transfer super-operator whose struc- 
ture is studied within a informational theoretical approach. Explicit expressions for computing the 
expectation values of symmetric observables are given both in the case of finite size systems and in 
the thermodynamic limit of infinitely many particles. 



I. INTRODUCTION 



The physics of strongly interacting many-body quan- 
tum system is central in many areas of physics. Our 
ability of simulating them is based on the possibility to 
find an efficient description of their ground state. This is 
the case, for example, of White's Density Matrix Renor- 
malization Group [H which can be recasted in terms of 
Matrix Product States (MPS) 0, i, i, i, 0. Such rep- 
resentations are characterized by a tensor decomposition 
of the many-body wave-function which allows one to ef- 
ficiently compute all the relevant observables of the sys- 
tem (e.g. energy, local observables, and correlation func- 
tions), and to reduce the effective number of parameters 
over which the numerical optimization needs to be per- 
formed. MPS fulfill these requirements and can be used 
to describe faithfully the ground states of not critical, 
short range one-dimensional many-body Hamiltonians. 
However MPS are not efficient in providing an accurate 
description in other relevant situations, i.e. when the 
system is critical, in higher physical dimensions or if the 
model possesses long-range couplings. Several proposals 
have been put forward to overcome this problem. Pro- 
jected Entangled Pair States 0] generalize MPS in di- 
mensions higher than one. Weighted graph states Q can 
deal with long-range correlations. Here we focus on a 
solution recently proposed by Vidal [9] who introduced a 
tensor structure based on the so called Multiscale Entan- 
glement Renormalization Ansatz (MERA). The MERA 
tensor network satisfies both the above efficiency require- 
ments and accommodates the scale invariance typical of 
critical systems [lo|, [HI • The relevance of this approach 
might represent a major breakthrough in our simula- 
tion capabilities [l2[ and motivates an intensive study 
of the MERA - e.g. see Refs. [H, [H, [H, [H, [T^. In 
Ref. some of us described a connection between the 
MERA and the theory of completely positive quantum 
maps [l8|. In the context of homogenous MERA's (i.e. 
MERA's formed by identical layers of tensors - see be- 
low for details) this permits to introduce a transfer ma- 
trix formalism in the same spirit as it has been done 
for MPS 0, [Bl, [l^l^ while providing new tools to compute 
physical observables using MERA's. As a result a connec- 



tion between the critical exponents governing the decay 
of two-points correlation functions and the eigenvalues 
of the MERA transfer matrix was identified yielding a 
simple method for determining the properties of critical 
many-body systems in the thermodynamic limit [1^, [13] • 



In the present paper we shall review some of the re- 
sults introduced in Ref. [l5[ providing explicit deriva- 
tions and clarifying the underlying mathematical aspects 
of the problem. In particular we formalize an important 
property of homogeneous MERA's by presenting two the- 
orems that allow one to evaluate the expectation values of 
symmetric observables (including translationally invari- 
ant Hamiltonians) in terms of a unique MERA transfer 
super-operator. Furthermore the thermodynamic limit 
of the MERA states is analyzed clarifying the condition 
under which such limit exists. 



The paper is organized as follows. Sec. [Ill is devoted 
to review the basics of the MERA tensor network. We 
discuss their causal cone structure and introduce the sub- 
set of homogenous MERA states. Even though most 
of this material can be found elsewhere 0, [2^ we de- 
cided to insert it here to make the paper self-consistent. 
This Section introduces also a new theoretical tool (the 
causal shadows of the MERA) which will play a fun- 
damental role in the subsequent derivation. Sec. [IIll is 
the central core of the paper: here we analyze the quan- 
tum channel description of MERA's showing how global 
quantities such us energy, average magnetization, etc. of 
a homogeneous MERA state can be described in terms 
of a single super-operator (the average QuMERA chan- 
nel). Also the MERA transfer operator is defined by 
moving in the Liouville representation [H, |2l| (the lat- 
ter is reviewed in Appendix [A]) . In Sec. HVl and Sec. [Vl 
we then discuss the thermodynamic limit of a MERA 
state and the scaling behavior of its two-point correlation 
functions by using general properties of mixing quantum 
channels [23, 24, 25, 26] and exploiting the spectral prop- 
erties of the associated QuMERA channel. The paper 
ends with the conclusions in Sec. [Vll 
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FIG. 1: a) Graphical representation of a typical one dimen- 
sional MERA tensor network for a many-body system of 
= 16 sites. The red elements correspond to the disentaglers 
tensors the blue elements are the isometry tensors A, and 
the green element C is the hat of the MERA. Any two joined 
legs from any two distinct nodes indicate saturation of the as- 
sociated indices of the corresponding tensor 0. The dashed 
line indicate periodic boundary conditions (i.e. the right-most 
X re-emerge on the left of the graph) . Alternative MERA de- 
compositions can be obtained by reordering the links of the 
graph, e.g. see Ref. [16]. b) Representation of the contraction 
rules of Eqs. (|2]), (|3]) imposed on x ctnd A (here the inverted 
elements represents their adjoint s counterparts). 



II. THE MERA NETWORK 

In this section we give a brief, self-consistent review of 
the basics of the MERA tensor decomposition which was 
introduced by Vidal in Ref. [9]. The only new element 
is the formalization of the notion the Causal Shadows 
presented in Sec. Ill C[ 



A. Basics 

Consider a many-body quantum system S composed 
by = 2^ sites of dimension d (qudits). Its pure states 
can be expressed as 

where for j G {1, • • • , N} and £ e {1, • ' ' ^d} the vectors 
G Hd form the computational basis of the j-th sys- 
tem site and where the type-( ^ ) tensor %i,£2,--- '•— 
(C^i 7 ' ' " 7 ^^AT 1^) cire the associated probability ampli- 
tudes. The MERA representation Q assumes a de- 
composition of T in terms of a collection of smaller, fi- 
nite size tensors, which differently from the linear MPS 
structure Q , are organized in a complex two-dimensional 
graph. An explicit example is shown in Fig. [H Here the 



links emerging from the lowest part of the graph repre- 
sent the N physical indices of T associated with the sites 
of S. The nodes of the graph instead represent tensors. 
They are divided in three groups: the type-( I ) disen- 
tangler tensors x of elements Xe^'e^ represented by the 
red Xs; the type-( I ) tensors A of elements A^^^^^ repre- 
sented by the blue inverted Ys; and the type-( ° ) tensor C 
of elements C^^, £2/3/4^ represented by the green blob. As 
shown in Fig.[T]the x's, the A's are coupled in together to 
form a triangular structure with C as the closing element 
of the top: any two joined legs from any two distinct 
nodes indicate saturation of the associated indices. Con- 
sequently, apart from the hat, T is written as a network 
of 0{N) smaller tensors organized in m = log2(A^) — 2 
different layers which we enumerate from the bottom of 
the graph. 

For a generic MERA the tensors entering the decom- 
position may differ from node to node and their indexes 
may have arbitrary (finite) dimensions. Under these con- 
ditions any state of 5* can be represented as in Fig. [1] by 
a proper choise of the x's and the A's. In the follow- 
ing however we will restrict the analysis to the special 
class of homogeneous MERA states in which all the x's 
and the A's entering the decomposition are identical and 
in which all the indexes of the graph have the same di- 
mension d [22^. With this choice the MERA identifies a 
much narrower but more treatable subset of many-body 
quantum states. The interest in such subset is motivated 
by the fact that homogeneous MERA's possess an intrin- 
sic scale invariance symmetry built in which is typical of 
critical, translationally invariant systems. In particular 
by removing the first layers from a A^-site homoge- 
nous MERA state |^) we obtain smaller versions of such 
vector constructed with only sites which, how- 

ever, in the limit of sufficiently large N (thermodynamic 
limit) still preserve the same correlations of the origi- 
nal one. Such a symmetry is believed [1, [H, [IB, [l^ to 
be sufficient for characterizing (at least approximatively) 
the ground state properties of critical, translationally in- 
variant Hamiltonians, making homogenous MERA states 
optimal candidates for their numerical simulations. 

B. Causal cones 

What really makes the MERA decomposition a conve- 
nient one is the assumption that the tensors composing 
the graph satisfy special contraction rules [9] - see Fig. [1] 
part b). Specifically one requires the following identities 

[x-x]TS = l^-^^TS = CC' (2) 

[A -A]; = (3) 

where S is the Kronecker delta, represents upper- lower 
contraction of consecutive tensors [27] , and where x 
A are the adjoint s of x and A defined by 

-ui,U2 _ (^txM >|* \ni,n2 _ /A^i \* (A\ 
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Expressed in operator language, Eqs. ([2]), (|3]) imply that 
X and A can be interpreted, respectively, as unitary trans- 
formation acting on two qudit sites and as an isometry 
that maps one qudit into two qudit s. 

Under these constraints each triple formed by three 
consecutive sites of the system is associated with a causal 
cone (CC) identified via percolation - see Fig. [2j Only 
the x's and the A's belonging to the CC can contribute 
not trivially in the evaluation of the of expectation val- 
ues of the local observables acting on such triple [1^. 
This is an important property of the MERA which al- 
lows one to reduce the number of contractions that need 
to be performed when evaluating expectation values on 
1^) from 0{N log2 N) to only 0(log2A^), exponentially 
simplify the complexity of the calculation Q . As a result 
given Aj^ an observable acting not trivially on the triple 
of sites k — l^k and /c + 1, we can express its expectation 
value as 

mAum = (4 • qI'")) • A • {Q^r^ ■ cu) , (5) 

where Ak is the tensor associated with the operator 
-^k [l^, is tensor associated with the CC of 

the triple /c, and Ck is the hat tensor with k specifying 
which of its 4 lower indexes couple to Q^^^ (as before 
and the ".7." represent upper-lower index contraction 
and the adjoint operation). Q^^^ is obtained by properly 
cascading m copies of the following type-(^ I ) tensor A^, 

where, as in Ref. [l5[, for easy of the notation we use 
typographic symbols o, o, ^ to indicate summation over 
the corresponding index. The way such tensors couple 
with each other and with Ak and Ck is specified by their 
position within the cone and ultimately depends upon 
the location on the triple in 6*. As shown in the figure 
there are two possibilities. Specifically the Ai tensor of 
the (m' — l)-th layer of the CC can saturate its upper 
indexes ui^U2^us either with the lower indexes -^2,^3,^4 
(modality L) or with the lower indexes ^3, ^4, ^5 (modality 
R) of the M. tensor of the m'-th layer. We distinguish 
the two cases by assigning a label a e {L^R} to each 

element of Q^^^ - see Fig. [2j)/ With this choice we can 
now write 

Qt^ = A^^<.) ■ M^,.> ..... , (7) 

where, for G {1, • • • , m} and /c G {1, • • • , N} the index 
a[^/ G {L^ R} specifies which lower indexes A4 (k) uses to 
connect with A4 (k) (or with Ck and Ak if = 1 or 

m,' — 1 

= m). 

Analogous simplifications occur also for non-local ob- 
servables. Of particular interest are the 2-point correla- 
tion functions of the form {"^{Ak Bk'l"^) with Ak and 
Bk' being (local) operators which act not trivially on the 
triples formed by the sites /c — 1, /c, k -\- 1 and /c' — 1, k' , 




FIG. 2: a) The black elements of the graph represent the ten- 
sor Qk of the causal cone associated with the triple formed by 
the 5-th, 6-th and 7-th sites of the system - see Eq. (|7|). It can 
be identify by connecting the triple with the MERA's hat via 
percolation starting from the bottom of the graph. Thanks 
to the contraction rules ([2]) and Q the remaining x's and 
A's (in gray) do not contribute when evaluating the expecta- 
tion values of observables which act locally on the triple. The 
light box underlines the tensor Ai of Eq. Here N = 16. 
b) The two alternative ways in which a tensor M can enter 
in Qk'. the empty circle represent the links that connect with 
the neighboring elements of the cone. 



/c' + 1, respectively. In this case the contraction rules ([2]), 
([3|) determine a joint CC for the sites /c, k^ formed by two 
single-triple CC (one for each triple), which intercept at 
the m + 1 MERA layer (counting from the bottom of the 
graph) with [l5[, 

m = int[log2 \k - k'\] - 1 , (8) 

see Fig. [3l This allows us to express (g) Bk'l"^) as 

{X,,, . (Qf ^Qif^))) . (AM . {(Qt^Qp^) . ,(9) 

where AkBk' is the tensor associated with Ak Bk' ^ 
Sfc^^Si^^ is tensor which describes the CC up to the 
m-th layer (it is given by the product of two indepen- 
dent single triple CC ([7j)), and where Xkk' describes the 
convolution of the MERA hat with the remaining part of 
the CC (i.e. the part above the m-th layer). 

C. Causal Shadows 

A notion which is complementary to CC is provided by 
what we define the Causal Shadows (CS's) of the MERA. 
If CC's play a fundamental role in the calculation of ex- 
pectation values on the MERA state, the CS's are funda- 
mental in simplifying the analysis of symmetric quantities 
as will be clear in the next section. 

Given a certain set of links C of the m'-th MERA lay- 
ers, we define its associated CS as the set of all CC's that 
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causal cone k causal cone 



FIG. 3: The black and green elements represent the causal 
cone structure associated to the product Ak Bj^f (here N = 
32, /c = 9, and k' = 24): it is formed by merging the CC of the 
k triple (black elements) with the CC of the k' triple (green 
elements) which intercept at the (m+l)-th MERA level. The 
empty circles describe the 6-qudits quantum state on which 
the two single-triple CC's operate upon. It is associated with 
the tensor Affefc/ of Eq. ^ and with the density matrix ak,k' 
ofEq. dSni). 
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FIG. 4: a) The black links represents an example of first 
neighboring causal shadow of depth 2 associated with three 
consecutive links of the second layer oi a N = 32- MERA. 
b ^Tensor of a causal shadow of depth 2 associated with a triple 
The insets show its decomposition in terms of products (|7|): 
starting from the top-left corner and moving clockwise we 
have Ml • Ml, Mr • Mr, Ml - Mr and Mr • Ml- 



allow one to reach elements of C (and only those) from 
the physical indexes of the MERA (i.e. the bottom of the 
graph) and define its depth - more precisely CS is the 
set of x's and A's belonging to such CC's. A trivial ex- 
ample of CS is obtained by considering m' = m (upmost 
MERA layer) and identifying C with the set of 4 emerg- 
ing links: in this case the CS includes all x's and A's of 
the MERA. Less trivial examples are shown in Figure [H 



For each CS we can clearly associate a tensor (this is the 
tensor formed by the x's and the A's of CS) and a set of 
physical indexes (this is the set of physical indexes that 
are attached to the CS tensor). Of particular interest for 
us are the CS's associated with triples of links as those 
shown in Fig. [H A simple combinatorial analysis shows 
that the physical indexes of such causal shadows contain 
2^+2 contiguous elements with m' being the depth of 
the CS. Most importantly one can verify each one of such 
CS tensor contain all possible 2^ sequences formed by 
combining m' tensors M either with modality L or with 
modality R: this is a trivial consequence of the fact that 
the CS includes all possible paths ([7j) which ends in the 
selected triple — see Fig. |4]b) for an example. 



III. QUMERA CHANNELS 

A better insight and more compact expressions for the 
expectation values on homogeneous MERA's can be ob- 
tained by moving in super-operator language 
Within this approach the tensors entering the causal 
cones associated with a specific subset of MERA's sites 
are re-organized to form concatenations of certain quan- 
tum channels (the QuMERA channels) whose definitions 
and properties depend explicitly upon the selected subset 
of sites. In the following we will review this approach and 
provide an explicit proof of two Theorems that allow one 
to derive a simple analytical expression for the average 
quantities computed on the state |^) associated with the 
MERA. 



A. Local observables 

Let us start considering the case of observables Ak op- 
erating on the triple k formed by the neighboring sites 
k — l^k and k -\- 1 for which Eq. ([5j) applies (single- 
sites observables are trivially included as a special case). 
We notice that, depending on the value of a, the Ma 
tensors contributing to ([7j) can be associated with two 
families of operators {Lr}r and {Rr}r acting on the 
Hilbert space Hf^ and labeled through the composed 
index r := (ri,r2,r3) with ri^2,3 being d-dimensional. In 
the computational basis they are defined by the matrices 

and are related through a reshuffling H of the input and 
output qudits, i.e. 

Rr = U{Lr) := PLrP"^ , (12) 

where P = P''' is the unitary transformation which ex- 
changes the first and the third qudit. Most importantly, 
according to the contraction rules ([2]), (|3|), the sets {Lr}r 
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and {Rr}r satisfy the following normalization conditions 
^L.Lt = i®' = ^i?,i?t ^ (13) 

r r 

with i being the identity operator of Hd- Therefore 
{Lr}r can be used to identify a completely positive, 
unital, not necessarily trace preserving super-operator 



^^H^ which transforms the linear operators A of 
according to the following expression 



03 



(14) 



Analogously {Rr}r defines the map <l>^^ which is related 
with through the identity 



H 



n o o n 



(15) 



indicates the composition of super- 
operators [3^. We will refer to such operators (and to 
their compositions) as QuMERA channels (i.e. quantum 
MERA channels). We also introduce the following vector 
of Hf\ 

l^hat) := ^ ^^1,^2,^3,^4 \^il1^i2 1^i3 1^u) 5 (16) 

^1,^2,^3,^4 

which can be assumed to be normalized (thanks to 
Eqs. ([2]), (|3|) the norm of |^hat) and |^) coin- 
cides). For j G {1,2,3,4} we then introduce : = 
Trj[|^hat)(^hat|] the reduced density matrices obtained 
by tracing |^hat)(^hat| over one of its j-th qudits, e.g. 



^1 ^2,^3,^4 ^2'^3'^4 
^\^£2^ ^^3^ ^^4) i^i'^ ^i'^ 5 



(17) 



(assuming |^hat) to be symmetric under permutations 
the p^^'' becomes identical). With these definitions one 
can finally cast the expectation value ([5]) as 



(\E'|ife|\E')=Tr[pg'=)$(j-^ 



...o$^"^'')(i)], (18) 



where jk G {1,2,3,4} and where (enumerating from the 
lower MERA level of Fig. [J) <l>^^ is either the map <l>^^ 

or <l>^^ associated with the corresponding element Ma of 
the causal cone ([7|). Notice also that in the rhs we have 
removed the label k from since there is no longer 
the need to specify over which triple the operator is act- 
ing (the trace runs in fact over 3 qudits only). The ex- 
plicit value of jk as well as the sequence of maps entering 
Eq. (p!8|) depend upon N and k: for instance in the case 
shown in Fig. [2] (i.e. N = 16, k = 6) the super-operator 

sequence is o while j = 4. Exploiting Hilbert- 
Smidth duality Eq. (p!8|) can also be written as 



TT[^i4'^) o . . . o )(pg'^^) A] , (19) 



with being the super-operator in Schrodinger 
picture (i.e. its adjoint with respect to the Hilbert- 
Smidth product). This is a completely positive, trace 
preserving (CPT) channel [I^ whose operator sum rep- 
resentation is provided by the operators {Ll}r (if a = L) 
or by the operators {Rl}r (if a = R). Since Eq. ([19]) 
holds for all observables A^ we can finally conclude that 



$(<')o...o$(-J^'')(pg'=)) 



pk 



(20) 



with pk being the reduced density matrix of |^) associ- 
ated with the triple k. 

Building up from these results we now present 
a theorem which formalize previous observa- 
tions [i5|,[ii,[i3,[23: 

Theorem 1: Let A'^^^ be the symmetric version of the 
local operator A, i.e. A^^^ := j^Ylk^i-^k- ^ts expecta- 
tion value on the homogeneous MERA state |^) can be 
computed as 



=Tr[<l>^(pc) A] 



(21) 



where <l>^ := ^ o ^ o • • • o $ with $ being the equally 
weighted mixture of 



(22) 



and PC := ELi PcV^- 



Proof: From the definition of A^^^ and from Eq. ([2Q|) 
we can write 



1 ^ 



1 ^ 

i^Tr[^(^i'')o...o$(«^'^)(pg^))i] . (23) 



TV 



We now notice that when varying /c, <l>^^i o • • • o <|)('^1^J^) 
spans all possible m-long sequences of <l>*^^^ and <l>(^). As 
a matter fact each of such sequence is counted four times 
(one for each possible values of jk)- Remembering that 
m = log2 TV — 2 we can thus reorder the summation over 
k on the last term of Eq. ([23|) as follows 



^ ^ ^ Tr[$(»^) o . . . o $('»'")(^g)) A] , (24) 

where a is the string (ai,a2,--- ,am). Using then the 
identity 



ae{L,R}^ 

this finally gives Eq. ([2T]) . ■ 
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Remark:- An alternative proof can be constructed by 
expressing the involved tensor contraction in terms of 
the CS's associated with all possible triples of hnks of 
the MERA hat and by exploiting the fact that each of 
such CS's contains all possible combination of J\4l and 
Mr (i.e. of and <l>(^)). 

An important application of this theorem is obtained 
by considering the expectation value of translationally 
invariant Hamiltonians H with first nearest- and second 
nearest-neighbors coupling, i.e. 



H := 



N 



A'" 



^ (S) ^ (2) 

with H!>_\ ■ i+i describing 3-body and 2-body in- 

teractions, and with Hi being local terms (here i = and 
i = + 1 are identified with i = N and i = 1 respectively 
to enforce the proper periodic conditions). This can be 
expressed as 



N 



(27) 



k=l 



where 



(3) 

/c-l,/c,/c+l 



(2) 

/e-l,/e 



(2) 

/c,/c+l 



(1) 

k-l 



(1) 



H 



(1) 

/c+l 



(28) 



is the Hamiltonian terms associated with the triple 
formed by the sites k — 1, /c, and /c + l. We can thus 
interpret H /N as the symmetric version of the local ob- 
servable h. Therefore from Eq. ([2T]) follows the identity 



{mm 

N 



Tr[$™(pc) h] . 



(29) 



This expression shows that, for homogeneous MERA's, 
the evaluation of the average energy per site E = 
(^|i^|^)/7V of a generic translational invariant Hamil- 
tonian H can be expressed in terms of the channel <l> 
and of the (symmetric) reduced density operator of the 
MERA's hat. 



B. Two-points correlation functions 

Let now focus on the correlation functions of the form 
{"^lAk^Bk']"^) with Ak and Bk' being generic observable 
operating on the k ^ k' triples respectively. Applying the 
derivation of the previous section to Eq. ([9]) we can write 



= Tr[$^^' 



(30) 



where for G {!,••• , m} the m' ^ are 6-qudits 
QuMERA channels (Schrodinger picture) associated with 
the first fh levels of the casual cone of the sites /c, k^ (see 
Fig. [3]), while akk' is a 6-site density matrix obtained by 
"evolving" the MERA hat with the remaining part of 
the causal cone and tracing out some of the links (which 
one is indicated by the indexes /c, k' and ultimately de- 
pends upon the interception between the two indepen- 
dent causal cones of k and k^ - see Fig. [3] for an example: 
here the 6-sites of akk' are indicated by empty circles). 

An explicit expression for m' Ms obtained as follows 



(k,k'). 



(31) 



with <|)^^rr.'^ 



and ^(^-') 



being the single triple maps as- 
sociated with the causal cones of k and k' respectively. 

/■i{k,k ) \ 

Thus, depending on /c, k' the map rn' ' will be one 
of the following four channels, ^^^^ (g) (^^^\ ^^^^^ $(^)^ 
<|)(^) (g) <|)(^) or <l>(^) (g) As in the case of Eqs. ([19]), 

([2Q|) we can then use the fact the Eq. (|3Q|) holds for all 
possible two-triple observables to conclude that the joint 
state of the triples /c, k' can be expressed as 



/,(fc,fc')x /u(k,k')^ 



Pk,k' 



(32) 



with pk^k' the reduced density matrix of the MERA state 
1^) associated with such triples. 

The evaluation of Eq. ([3Q|) is in general quite compli- 
cated as it requires to compose four different maps in a 
specific order determined by the involved CC's (of course 
for some clever choice of k and k' such a sequence could 
be relatively simple to compute). One would be tempted 
to solve this problem by "symmetrizing" the two-point 
operator as in the local observable case (e.g. replacing 

Ak Bk' with Y^k^i Ak Bk^Ak , where Ak = k' -k). 
Unfortunately this is not sufficient, the reason being ul- 
timately related with the fact that even for homogenous 
MERA the state |^) is in general NOT translational in- 
variant for finite [311]. 

One way to circumvent this is to exploit the CS struc- 
ture to enforce a "local" symmetrization of Ak Bk'- 
Specifically, let A being a generic observable acting on 
a triple of qudits. Consider then a causal shadow CS^ 
of depth fh characterized by M = 2^ + 2 physical in- 
dexes {4^,4^+1, ••• ,4^+^}_ which win be grouped 
in a sequence of consecutive M — 2 triples labelled as 
/ca + 1, /^^A + 2, • • • , /ca + M — 1 (here kA is the leftmost 
physical index of CS^ whose explicit value depends upon 
the position of the CSa within the MERA — see Sec. Ill CI 
and Fig. [5]). We define the shadow operator as the aver- 
age of A over the triples of CS^, i.e. 



/CA + 1 + 2" 

E 

k=kA-\-l 



Ak 



(33) 



where Ak is the operator A acting on the k-th triple of 
CSa' Consider then a second operator B and a second 
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causal 
shadow CSa 




causal 
shadow CSb 




indexes of the two causal shadows - see Fig. O 

Proof: The proof proceeds as in the case of Theorem 
1. First of all we use Eq. (|36|) to write 



^Tilpkk' {Ak^Bk^)] ,(38) 



FIG. 5: The black and green elements represent the causal 
shadows CSa and CSb respectively. The arrows indicate their 
leftmost sites kA and ks (here N = 32, kA = 8, and ks = 20). 
The empty circles describe the 6-sites described by the density 
matrix a^'^^'^'^^ of Eq. (l37l) . 



CSb of depth m which is first neighbor with CS^ (that 
is the rightmost index at the top of CS^ is first neighbor 
with the leftmost index at the top of CS^, or vice- versa 
- see Fig. |4]). It is worth noticing that the physical sites 
of CSa and CS^ are separated by a distance which is 
exponentially large in m. Indeed such distance can be 
easily computed as 



kA = 2{2^-l). 



(34) 



with ks being the leftmost physical triple of CS^ (here 
for the sake of simplicity we assume ks > kA)- Another 
important property is the fact that given CCa and CC^ 
causal cones associated with CSa and CS^, respectively, 
they are independent up to the level fh but intercept at 
the level m + 1. As a matter of fact one can easily verify 
that all couple of CC's that share this property enters in 
CSa and CSb- 

Also for CSb we define now a shadow operator 



Kb 



1 



and consider the joint shadow observable 



^kA 



(35) 



(36) 



This is a symmetrized version of Ak Bk' obtained by 
averaging locally over all possible choices of k and k' 
whose causal cone intercept in the same points of the 
MERA. 

Theorem 2: Let AB^^ be defined as in Eq. (|36|) . Its 

expectation value on the homogeneous MERA state |^) 
can be computed as 



(^I^^L,Aj^)=Tr[(^-0^-)(a 



(fcA,Afc) 



)(i®B)p7) 



where fh is the depth of the causal shadows CSa ctnd 
CSb, ^ cis in Eq. ([22]) . and where now ai^^Ak) 
the 6-qudit density operator associated with the upper 



where the summation is performed over the physical sites 
of CS^ and CSb and where pk^k' is the joint reduced den- 
sity matrix of |^) associated with the triples k and /c^ 

.,(fc,fc')\ 

The latter can be expressed as in Eq. (|32|) with m' ^ 
being the 2-site QuMERA channel (|3T]) associated with 
the causal cones of the triples k and k\ Most importantly 
in our case the 6-qudits density matrix ak^k' is indepen- 
dent from k and k^: in fact it represents the state of the 
(m + l)-th level that is attached with the CC's associ- 
ated with the triples k and k' and by construction all 
the CC's belonging to a given CS intercept in the same 
points. Exploiting this we identify ak^k' with a^^^^^k) Qf 
Eq. (|37|) . Furthermore the (independent) average over 
k and k^ allows one to generate all possible m-long se- 
quences of <l>^^^ in both CS's. We can thus write 
the rhs of Eq. ([38l) as 



k,k' 



(2^)2 

x((j('=-*''^'=))(i(g)B)] 

— -"-^[/^ 2^ ^ 2^ 2^ 

k k' 

x(a(^^'^^))(i0^)] . 

Equation (|3Q|) finally follows by the identity ([25|) . ■ 

Remark:- We can further symmetrize the correlation 
function by averaging Eq. (|36|) with respect to their ab- 
solute position within the MERA (keeping the relative 
distance among the CS's constant), i.e. by replacing 

^^Ha, with 



1 ^ 



Afc ' 



(39) 



kA = l 



where kA is the leftmost triple of CSa- With this choice 
Eq. ([371) still applies by substituting a^^AAk) ^^^h its 
average counterpart a := X^^^^ (j^^^'^^^A/", i.e. 



(^I^^aII^) =Tr[(^^0^^)(a) {A^B)] . (40) 

In Sec. [Vl we shall see how Eq. (|^Q]) can be used to 
determine the scaling behavior of two-point correlations 
function of a MERA. 
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C. The MERA transfer operator 

Theorem 1 and 2 formany show that in extracting local 
or correlated (average) quantities of |^) one can focus 
on a single CPT map ^ obtained by averaging over all 
possible triple-sites QuMERA channels. This leads us to 
identification of a transfer operator for the MERA [15] in 
close similarity to what happens in the case of MPS (e.g. 
see Ref. [19^). 

The idea is to move to the Liouville space representa- 
tion [H, [2l| in which by "doubling" the degree of freedom 
of the system, the super-operators are represented by ma- 
trices, and the operators by vectors (see Appendix [Al for 
details). Specifically this is done by fixing an orthonor- 
mal basis {|^); ^} [32] on the three qubits space Hf^ and 
associating to each operator A of such system a vector of 
|i)) enf^Hf defined by 



\A)) :=^(z|i|0 



(41) 



According to this formalism the QuMERA channel <l> of 
Eq. ([22l) can now be described by the transfer opera- 
tor [liracting on ^ ^ defined by 



(42) 



where we used Eq. (|A6P of the Appendix and the fact that 
{L\/ \/2^ RI/ \/2 ; r} is a Kraus set for ^. Consequently 
we can exploit the identity ()A8p to write Eq. ([2T]) as 



(43) 



Similarly we can proceed for the two-point correlation 
functions (|iQ|) by defining the transfer operator of the 
channel acting on two-triples. Constructing the Li- 

ouville space as the tensor product of the Liouville spaces 
of the two triples the latter can be expressed as E^ , 
while Eq. (j30|) becomes 



{{A^B\{E^r®{E^r\G)) . (44) 



IV. THERMODYNAMICAL LIMIT 

In this section we analyze the property of homoge- 
nous MERA states in the thermodynamical limit of N ^ 
oo. To approach this problem we introduce the family 
^(x,A,C) composed by MERA states of exponentially 
increasing size, i.e. 

^(X, A,C) := {\^n) : iV = 2^ forn ^ 3 integer} , (45) 

where I^at) are MERA states with N sites constructed 
with the same tensors A and C. For each one of such 
family we can then use the result of the previous section 



to compute the thermodynamical limit of the expecta- 
tion values of (symmetrized) local observables as follows 



^(th)(^) := lim (^Ar|i(^)|^Ar) 

N^oo 

= lim Tr[$™(/5c) A] , 



(46) 



with A^^^ and <l> as in Sec. Illll and where in the last term 
we used the fact that for a A^-sites MERA m = log2 N—2. 
The existence of 74^*^^(^) depends thus on the conver- 
gency of the limit lim^^oo ^^(pc), with <l> being the 
average QuMERA channel defined by the tensors x and A 
of the selected family ^. This allows us to relate the ther- 
modynamical limit of homogeneous MERA's with the 
problem of the convergency of repeated applications of 
a given CPT map and thus ultimately with its mixing 
(or relaxing) property [H, [H, [H, [2^ . 

It is worth reminding a mixing channel <l> is character- 
ized the property 



lim 

m^oo 



(47) 



with being the CPT map which (times Tr[5]) trans- 
fers every operator B into a fix density matrix pf (the 
fix point of <l>). In the following we will indicate ^/ as 
the final point channel of <I>. It satisfies the following 
important property: 



{^f^lY){&XY) =Pf^&Y 



(48) 



where Qxy is a joint operator of the system X on which 
<!>/ upon (i.e. three qudits) and of a generic ancillary 
system Y; Xy is the identity map on Y] and finally Oy : = 
Trx[Oxy]- It is a well know fact that the vast majority 
of CPT maps acting on given system are mixing (the 
non- mixing one form a subset of zero- measure). This 
clarifies that, a part from some rare pathological case, 
the limit (|46|) is well defined [[33] . Furthermore it gives us 
a simple way of computing such quantity. Indeed owing 
to Eq. (|47|) we can write 



A(*^)(^) = Tr[p/ A] 



(49) 



with pf being the fix point of the QuMERA channel <I>. 
The latter can be easily determined by solving the eigen- 
value problem 



Hp) = p , 



(50) 



as for mi xing maps is the unique solution of such an 
equation [IJlli, |25|, |2^ . 

More generally we have the following statement 

Lemma: The fix point pf of the QuMERA channel 
^ of the family ^(x,A,C) (when defined) coincides with 
the thermodynamical limit of the average reduced density 
matrix of the triple of the associated system, i. e. 



N 



Pf 



lim Vpfc/A^ , 

N^oo ^ — ' 



(51) 
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with pk being the density operator associated with k-th 
triple 0/ 1^). 

Proof: exploit the fact that Eq. (|49|) holds for all 
observables A and the fact that the expectation value of 
can be expressed as an average over all triples of 
the system — e.g. see the first line of Eq. 



Remark: Since <l> only depends upon x A, the 
average reduced density matrix of the family (as well as 
the quantities (|49|) ) does not depend upon the MERA 
hat tensor C. 

In a similar way we can also compute the thermo- 
dynamic limit of two-point correlation functions ([39|) in 
which we keep the distance Aj^ constant. In particular 

given AB^ as in Eq. (|39|) we define 



lim (*yv|AB^J*jv) 



(52) 



= lim Tr[($™ (g) $™)(a) {A ® B)] 

N^OO 

= Tr[($™ ® $™)(ct(*'')) (i ® B)] , 



with 



a^*"^) = lim (T , 

N^oo 



(53) 



being the thermodynamical limit of a, the latter being 
defined as the average reduce density matrix of 6 con- 
secutive sites associated with (m + l)-th MERA layer. 
Notice that differently from ([46|) the limit N ^ 00 does 
not translate into an infinite sequence of applications of 
the QuMERA channel ^. This is because the latter de- 
pends only upon the interception among the CS's inter- 
ception which is fixed by the distance A/^. On the con- 
trary varying N implies a variation on a which is taken 
into account by Eq. ([53|) . By exploiting the scale invari- 
ance of homogenous MERA's one can easily verify that 
a coincides with the thermodynamic limit of the average 
reduced density matrix of 6 consecutive physical sites of 
the MERA. Therefore it is possible to show that explicit 
expressions for a*^*^^ can be obtained by applying proper 
(multi-site) QuMERA channels to pf. In the following 
section we will not discuss this topic any further, instead 
we shall focus on the scaling behavior of Eq. ([52]) in the 
limit of large distances A^. 



V. SCALING BEHAVIOR OF TWO-POINT 
CORRELATIONS FUNCTIONS 

The scaling behavior of the two-points correlations 
functions for a homogenous MERA can be determined 
by looking at the spectral properties of its QuMERA 
channel ^ [H, [H, [l^. This can be done both for fi- 
nite dimensional MERA's and in the thermodynamical 
limit thanks to Eqs. (|4Q|) and ([52|) which recast the com- 
putation of the correlation functions in term of similar 



expressions. Here we will focus on the latter case which 
is by far the more relevant. 

First of all, given A, B generic observables acting on 
triple of sites, we introduce the following rescaled quan- 
tity 

AAB^l^\^) := AB^l^\^) - A^'^\^)B^'^\^) (54) 
= Tr[(<l>^ <l>^)(a(*^)) (Ai A^)] , 

where we used the fact that the 3-site reduced density 
matrix of a^^^^ is pf and where AA := A - A'^*^)(^), 
A^ := 5-5(*^)(^). The idea is to invoke once more the 
mixing properties of <l> which guarantee that this operator 
has a unique unitary eigenvalue [l^, UBI, • Exploit- 
ing then the spectral decomposition of ^ and keeping 
the highest order contributions this can now be written 
as [15i] . 



m>l 



c \r]r] 



log2 17777'! 



(55) 



with 77, 77' 7^ 1 being the eigenvalues of ^ of largest mod- 
ulus which contribute non trivially in the expansion, and 
where in the last term we used the fact that fh scales 
logarithmically with the distance A^ as in Eq. In 
this expression c is a term which scales at most polyno- 
mially on tti, i.e. c ^ 0(Poly{m)) :^ (9(Poly(log2 A/.)) 
— see below. Equation ([55|) shows a polynomial decay of 
the two-point correlation function of the system which is 
typical of critical system [35[. Its derivation resembles a 
similar calculation performed in Ref. [l^ for MPS's. It 
can be obtained by expressing AAB^^^ (^) in the Liou- 
ville representation which as seen in Sec. IIII CI gives 



AAB^ll\^) = {{AA^AB\{E^Y 



.Hf)^' which rep- 



where \a^^^^)) is the vector of {Uf^ ^ 
resent the state This can now be simplified by 

means of the identity Eq. (|A2ip and observing that for 

each vector \a)) G (Hf^ ^Hf^)^^ the following identi- 
ties applies 

{{AA^AB\{Ef^i)\a)) = {{AA0AB\{^^I){a))) 
= {{AA^AB\pf0p)) = {{AA\pf)){{AB\p)) 
= Tt[AA Pf] Tt[AB p] = , (57) 

(here p is the reduced density matrix of a and we used the 
property Eq. (gll) of <l>, and the fact that Tt[AA pf] = 0). 
Similarly one has 

((Ai A^|(i Ef)\a)) = . (58) 

Exploiting these identities we can now write 

A^Bi*|;)(*) = ((Ai® AB|(A^^® A£:™)|a("^))) 

■ (59) 
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where the being trigonometric, polynomial func- 

tions of fh denned by 

Cp := {{AA AB\ej{m) e,v(m)|a(^^))) , (60) 

with ejj'{fh) as in Eq. (|A23p . Equation ([55|) finally fol- 
lows by taking the couple j, / which has the largest value 
of \r]jr]j^\ and for which cjj^^ ^ 0. 

A. Self-adjoint transfer super-operator 

Of special interest is the case of MERA's which have a 
self-adjont $ QuMERA channel ([22]) (i.e. ^ = with 
representing $ in Heisenberg picture). An example 
of such MERA's has been recently studied in Ref. [iQ] 
in the calculation of the ground state properties of Ising 
and Pootz model. In this case we can write 

$(i) = ^r?,Tr[e]i]e,- (61) 

j 

with r]j being the (real) eigenvalues of ^ and &j being 
the corresponding eigen-operator properly orthonormal- 
ized with respect to the Hilbert-Smidth scalar product (in 
particular if <l> is mixing then 770 = 1 is non-degenerate 
and Oj=o = Pf)- Therefore Eq. ([56]) yields 

AAB^ll\^) = ^ vfvjf Tr[(6]0 6],)(^^*'^)] 

X Tr[(6^- (g) e^v)(Ai (g) A^)] , (62) 

where the properties ([57|) and ([58]) has been used to re- 
move the contributions in j = or = from the 
sum. The above expression coincides with Eq. ([59]) 
by identifying cj^^ with the coefficients Tr[(6] 

et,)(a-(th))]Tr[(ej Of){AA A^)] times a phase fac- 
tor. A direct proof of this can easily be obtained by 
observing that in this case the Liouville representation of 
^ is provided by the Hermitian operator 

E^ = J2v,\0,)){{e,\, (63) 

3 

with being orthonormal. Under this condition the 

power-law scaling ([55]) becomes exact as the coefficient c 
is now independent from fh. 

VI. CONCLUSIONS 

In this paper we have presented an extensive analysis 
of the properties of homogeneous MERA states based on 
the quantum channels approach introduced in Ref. [l5[. 
In particular we have proved some Theorems which al- 
lows us to characterized the (local) average properties of 
such states in terms of the spectrum of a single QuMERA 



channel ^: Theorem 1 establishes that the expectation 
value of any (average) local observable can be evaluated 
on the fix point (eigenvector associated with the maxi- 
mum eigenvalue) of <I>; Theorem 2 instead gives an ex- 
plicit expression for the two-points correlation functions 
of the system. Both Theorems holds also in the thermo- 
dynamic limit of MERA states associated with an infinite 
number of sites. In particular Theorem 2 allows one to 
identify the (power law) scaling behavior of the MERA. 
For the sake of simplicity the analysis has been performed 
assuming a specific MERA decomposition but it can be 
trivially generalized to any possible variation of the lat- 
ter. 
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APPENDIX A: LIOUVILLE SPACE 
REPRESENTATION 

As anticipated in Sec. IIII CI the Liouvillle representa- 
tion is constructed by fixing an orthonormal basis 
on the Hilbert space of interest H (e.g. the three qubits 
space Hf^) and by defining the following linear map- 
ping from the space B{H) of the linear operators of H to 

A = - \A)) := E(^I^IOI^) ^ lO(Al) 

Simple but useful properties of the mapping (|Aip are the 
following rules, 

\AbC)) = {A^C^)\B)) , (A2) 
Tr[it^] = {{A\B)) , (A3) 

which hold for all operator A^B^C G B{l-t). In these 
expressions stands for transposition with respect to 
the selected basis while A'^ is the adjont of A - 

notice that one has {{A\ = [\A))]^ = T.^^ {i\A\i'y {i\^ {i'\ 
with {i\A\i')^ being the complex conjugate of {i\A\i'). 
It is also worth noticing that, according to the above 
expressions, the vector = \i) \i) satisfies the 
following identities 

Tr[B] = {{i\B)) , (A4) 
\B)) = (B®i)|i)) = (i®B^)|i)) . (A5) 

In this language a CPT map <l> operating on B{H) is 
described by an operator acting on H defined by 

:=^M, , (A6) 
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where Mg are a set of Kraus operators of ^. Equa- 
tion ()A6p is a consequence of (|A2p by noticing that for 
all B eB{H) one has 

|^(^))) = Y.\Msmi)) = Y,Ms^M:\B)) 

s s 

= E^\B)) . (A7) 

This expression shows that ^ and its associated matrix 
have the same spectrum (i.e. the same eigenvalues). 
It also allows us to express the expectation values on 
evolved operators as matrix elements of as indicated 
by the following expression 

Tr[it^(^)] = {{A\E^\B)) . (A8) 

Finally it allows one to compute the successive applica- 
tion of a CPT map as follows: 

= (^^nS)) . (A9) 

An interesting problem is to determine the limit for 
m ^ oo of {E^)'^ when <l> is mixing (as in the case 
of the QuMERA channel case). From Eq. (|T7|) we know 
that this must be the transfer matrix E^j^ of the chan- 
nel which maps every operator into the fix point pf. 
According to the above definitions this implies 

E^,\A)) = T^lM \pf)) = im)) 1/5/)) . (AlO) 

Since this must be true for all vectors we can con- 
clude 

]\m {E^r = E^^ = \pf)){{i\, (All) 

m^oo 

which is consistent with Eq. (|49|). 

1. Some facts about the spectrum of $ and 

Consider the set S{^) := {rjj^j} of the eigenvalues 
of They are defined by eigenvector equations of the 
following form 

m,) = vA ^ E^lSj)) = vj\ej)) , (A12) 

where are the vectors associated with the eigenvec- 

tors Qj of ^ (for each rjj there can be more than one 
|0j))). Since the matrix E^ is generally not Hermintian 
the kj will be not real. However since <l> is CPT one can 
show that r]j belongs to unit circle [18] (i.e. \r]j \ ^ 1), and 
that 7^0 •= 1 is always an element of the spectrum, i.e. 
1 G S{^). Furthermore one has that since <l> is trace pre- 
serving then the eigenvectors associated to eigenvalues 
7^ 1 are traceless operators, i.e. Tr[6j] = (in particular 
they cannot be density matrices). Finally one can verify 
that if rjj G S{^) than also its c.c. is an eigenvalue of 
^, i.e. 77* G S{^) (too see this just take the adjoint of 
the left hand side equation of ()A12p and use the fact that 



^(e)t = ^(et)). More generally one can verify that the 
couple r]j and 77* have the same Jordan structure, i.e. 
their corresponding Jordan blocks (see below) will have 
the same dimensions. 

As already mentioned the spectral properties of ^ de- 
termine uniquely the mixing property of the map (2^, 
[25|,[2^. In particular, it is known that the map is ergodic 
(i.e. it has a unique fix point that satisfies Eq. (|^) if and 
only if r]o is non-degenerate (i.e. it has a unique eigen- 
vector). Furthermore, it is known that <l> is mixing if and 
only if 770 = 1 is the only eigenvalue with unitary mod- 
ulus and it is (non-degenerate), i.e. \r]j\ = 1 iff r]j = 1. 
We finally remind that ergodic maps are not necessarily 
mixing even though any mixing channel is necessarily er- 
godic, and that mixing channels are dense in the set of 
the CPT maps. 

2. Jordan block decomposition of E^ 

We have seen that in general E^ is not Hermitian: as a 
matter of fact, typically it will not be even orthogonal i.e. 
diagonalizable - see however the discussion of Sec. IV Al 
and Ref. [iHl • We can however still put it in Jordan form 
by similarity transformation, i.e. 

E^ = fjf-^ , (A13) 

where T is an invertible operator and where J is the 
Jordan form associated with E^. A part from a triv- 
ial permutation of the blocks, the operator J is uniquely 
determined as J := (BjJdj with Jd^ {r]j) being the Jor- 
dan block of dimension dj associated with the eigenvalue 
r]j of It is worth reminding that each eigenvalue r]j 
can have more than a single block: the total number of 
such blocks corresponds to the so called geometric mul- 
tiplicity of 7^j, i.e. to the number of linearly independent 
eigenvectors of r]j (i.e. to the dimension of the associated 
eigenspace). Finally the sum of the dimension dj of all 
the blocks associated with a given eigenvalue r]j is to the 
so called algebraic multiplicity of 7/j, i.e. the number of 
zeros of the characteristic polynomial associated with the 
solution X — r]j [3^] . For diagonalizable matrices one has 
dj = 1 and the geometric multiplicity coincides with the 
algebraic one. We remind also that the operator JdjiVj) 
can be written as the following matrix 

Jd,{Vj) = vAd,^Nd^ , (A14) 

with tdj being the dj x dj identity matrix and where Nd^ 
is nilpotent matrix which satisfy the condition {Ndj )^ = 
for all p ^ dj (specifically it is either a dj x dj matrix of 
I's above the diagonal or is the null matrix). 

It is interesting to observe that the matrix E^ satisfies 
the following condition 

K = Y1 ^Ms = S{E^)S'^ (A15) 
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where S is the swap operator which sends any opera- 
tor of the form A ^ B into B ^ A. This is unitary and 
Hermitian (i.e. S = S"^ = S'^). Therefore and 
are connected through a similarity transformation (i.e. 
they are mapped into each other by an invertible oper- 
ator): consequently E^ and E^ admit the same Jordan 
form decomposition. This is an important fact that tell 
us that given an complex conjugate couple r]j and r]j of 

eigenvalues of E^ , their associated Jordan block will have 
the same structure (i.e. dimension and possible degenera- 
cies). From this point of view hence E^ share some prop- 
erties the real matrix. In particular one can decompose 
E^ in the so called real Jordan form where the (complex) 
blocks associated with each couple of complex conjugate 
eigenvalues can be grouped together to form real " super- 
blocks" [U- Finally it is worth remembering that E^ 
is connected through similar transformation also with its 
transpose E^ (this is a general property of all complex 
matrices [sj), i.e. E^ = R E^ with R invertible 
but not necessarily unitary. Furthermore since E^ is con- 
nected through E^ by swaps operation (the derivation is 
as in Eq. (|A15P ) it follows that these all these matrices 
are similarly equivalent i.e. 

E^^EI^^eI^eI, (A16) 

(where stands for the similarity equivalence), and 
will have the same spectra (3Q| . 

The above expressions are extremely useful when com- 
puting successive application of <l> - see Eq. ()A9p . Indeed 
from Eqs. (IA13P and ()A14p one gets 

{E^r = fj^f-' =f(^ej [JdAVj)r)T-^ , (A17) 
with [Jd,{r]j)r = Q^'^\m) and 

Q^^H^.):= E ("T) ^7M^d,r, (A18) 



being bounded operators which are polynomial in m. 
This implies that for \rij\ < 1 one has 

lim [4(r,,)]™ = 0,,. , (A19) 

m^co 

with 0dj being the dj x dj null matrix. Therefore for <l> 
mixing we can write 

lim [E^r = t(Ji{Vo) ®i/o 0d,)f-' = ^$,,(A20) 

where we use Eq. (|Alip and the fact that for a mixing 
channel 7^0 = 1 is not degenerate and thus its correspond- 
ing Jordan block Ji{r]o) is the 1x1 matrix formed by the 
single element 1. Thus Eq. (|A17p can now be written as 
follows 

{E^r =E^^^^AmE^ , (A21) 
where A YYiE^ is a contribution that nullifies for m ^ 00. 
It can be expressed as 

A^^^ :=^|7/,re,(m) , (A22) 

with the matrix ej{m) being a trigonometric, polynomial 
function defined by 

e,(m) := f[e,yoe'"' "'^^^^'^ Q^^^X^W'^ • (A23) 

It is finally worth mentioning that for Hermitian E^ the 
whole analysis simplify. In this case in fact T is a unitary 
transformation, the T^j's are real, while Q^'^\r]j) = h • 
Under this condition ej (m) become independent from m 
and coincides with projector on j-eigenspace of E^. 
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